CMB Constraints on a Stochastic Background of Primordial Magnetic Fields 
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We constrain a stochastic background of primordial magnetic field (PMF) by its contribution 
to the angular power spectrum of cosmic microwave background anisotropics. We parametrize 
such stochastic background by a power-law spectrum with index ub and by its Gaussian smoothed 
amplitude Bx on a comoving length A. We give an approximation for the spectra of the relevant 
correlators of the energy-momentum of the stochastic background of PMF for any ub- By using 
the WMAP 7 year data in combination with ACBAR, BICEP and QUAD we obtain the constraint 
SiMpc < 5.0 nG at 95% confidence level for a stochastic background of non-helical PMF. We discuss 
the relative importance of the scalar and vector contribution to CMB anisotropics in obtaining these 
constraints. We then forecast Planck capabilities in constraining BiMpc- 



I. INTRODUCTION 

The origin of the large scale magnetic fields observed in 
galaxies and clusters of galaxies is an open issue of great 
importance in modern astrophysics (see 1] for a review). 
Primordial magnetic fields (PMF) generated in the early 
Universe could have been the seeds for large scale mag- 
netic fields and have left an imprint in the anisotropy 
pattern of the cosmic microwave background (CMB). A 
primordial hypothesis for generating the seeds amplified 
afterwards by adiabatic compression and dynamo - can- 
not be discarded [1| , also in light of recent observatiqns^of 
strong magnetic fields in galaxies at high redshift [^ [SJ . 

PMF with a comoving amplitude of several nG can 
leave interesting imprints on CMB anisotropics. A 
stochastic background of PMF is modelled as a fully inho- 
mogeneous component and its energy momentum tensor 
(EMT) - quadratic in the magnetic fields - is considered 
at the same footing as linear inhomogeneities in the other 
components and linear metric fluctuations. A stochas- 
tic background of PMF generates independent modes for 
all kinds of linear perturbations: there has been sev- 
eral studies for scalar |4l-[ll|. vector [l 4 -|l4 | and tensor 
[l3l llSl [la | perturbations. See Refs. [^, [l3| for studies 
which take into account all the types of contributions. A 
stochastic background of PMF affects also the statistics 
of CMB anisotropics, and not only its power spectra: be- 
ing quadratic in the magnetic field amplitude, the EMT 
is non-Gaussian distributed [I3 and therefore the bis- 
pectrum of CMB anisotropics can also be a useful probe 
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In our previous works [^, Q we have refined the com- 
putation of CMB anisotropics in presence of a stochastic 
background of PMF: we have computed the initial condi- 
tions for cosmological perturbations in the radiation era 
keeping into account only relativistic degrees of freedom 
and the correlators for the Fourier transforms of the EMT 
in presence of a sharp cut-off which mimics the damping 
scale due to viscosity [^, Q for few values of the spectral 
index ub of the stochastic background. 

In this paper we use and extend the above theoretical 
description to derive the CMB constraints on a stochas- 
tic background of PMF which can be obtained by current 
and future data. We therefore use a modified version of 
CosmoMC [20] connected with a modified version of CAMB 
2l| containing all the above features to constrain B\ - the 
amplitude of the magnetic field smoothed over a comov- 
ing scale A - and tib with the most recent compilation of 
CMB anisotropics data, therefore updating previous in- 
vestigations d,!!!,!!!]. We also forecast the Planck ^^ 
capabilities in constraining such a background of PMF. 

Our paper is organized as follows. In Section II we 
present the theoretical set-up needed for deriving the 
CMB constraints. As new results we give an approxima- 
tion for the PMF energy-momentum valid for any ns, 
we discuss the correlator between the energy density and 
the Lorentz force and we extend the regular initial con- 
ditions for cosmological perturbations in the relativistic 
regime in presence of a stochastic background of PMF 
9] including the correction due to matter, collecting the 
details in Appendix A, B and C, respectively. In Sec- 
tion III we discuss the constraints from WMAP 7 years 
data [Ullel, ACBAR 27], BICEP [H and QUAD [H 
on a flat ACDM model plus a stochastic background of 



PMF. We present the Planck capabilities in constraining 
B\ and ub in Section IV and we summarize our results 
in Section V. 



II. STOCHASTIC BACKGROUND OF PMF 

A stochastic background of PMF acts as a fully inho- 
mogeneous source to the Einstein equations. As usual, 
we assume the infinite conductivity limit in which the 
PMF simply scales with time as B{x,t) — B{x.)/a{T)'^ 
where a{T) is the scale factor normalized to ao = 1 to- 
day and r is the conformal time. We model PMF with 
a power-law spectrum Psik) ~ Ak'^'^. The two point 
correlation function for inhomogeneous fields is: 



(B,(k)i?*(k')) = (27r)35(k-k')(<5, 



kk)^ (1) 



where ub > —3 to avoid infrared singularities. Magnetic 
perturbations survive the Silk damping but on smaller 
scales PMF are damped by radiation viscosity. We model 
this damping introducing a sharp cut off in the spectrum 
at a damping scale ko U^i] which is much smaller than 
the Silk scale. To parametrize PMF amplitude, we choose 
to use the most common convention which smooths the 
magnetic fields with a Gaussian filter over a comoving 
scale A: 



Bl = 
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(2) 
As damping scale we adopt the proposal of Ref. |12j, in 
which k]j is function of Bx and the spectral index ub '■ 
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(3) 

The EMT of PMF is a source term for the Einstein- 
Boltzmann system and, being these quadratic in the 
Bj(k), the Fourier transfoms of the EMT components 
are convolutions Il3l. [lal : 
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^^dpPB(p)PB(|k-p|)(l 
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\Il(^Hk)\^ = ^JjpPB{p)PB{\k-p\ 

X (1 + 272 + 72/32), 
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(5) 



(6) 



where ji = p • (k - p)/|k - p|, 7 = k ■ p, 
/3 — k-(k — p)/|k — p| and fl denotes the volume 
with p < ko- 

The analytical exact results for Eqs. (4-6) are given 
for specific values of ub in our previous papers [8|, |9|. 



However, the expressions for generic Hb are rather com- 
plicated and cannot be used in a numerical implemen- 
tation. Therefore we fitted the analytical results with 
easier expressions to be inserted as spectra for the PMF 
in modified version of CAMB ^] and CosmoMC 2Q]. Since 
the spectral shape varies with the spectral index we had 
to divide the spectra in different index ranges. The first 
natural split is between indices greater and smaller than 
Ub = —3/2: this division is very natural since is required 
by the change in the infrared behaviour between the two 
ranges. In order to achieve the best accuracy with the 
fits and considered the wide range of spectral indices that 
we wanted to explore we decided to do a further splitting 
in the ub > —3/2 range between positive and negative 
spectral indices. In the end we result into three differ- 
ent spectral fits for each EMT component. In Fig IT] we 
show respectively for scalar energy density, scalar Lorentz 
force and vector anisotropic stress the results of our fits 
compared with the analytical results. We note how the 
fits are in excellent agreement with the analytical results 
and we refer the reader to Appendix A for the details. 
In the next two subsections we address other two issues 
related to the theoretical characterization of a stochastic 
background of PMF, recently addressed in Refs. |l(]ll30|. 



A. Cross-Correlation between Lorentz Force and 
Energy Density 

The first issue is represented by the cross correlation 
between Lorentz force and magnetic energy density as 
pointed out in [1^. In particular the analysis of the 
magnetized scalar mode involves three quantities: mag- 
netic energy density, Lorentz force and anisotropic stress. 
These quantities are not independent one from the other 
but are related throught the conservation equation for 
the PMF EMT (in the magnetic hydrodynamic limit): 
aB = ^ + Lb- We approximated the Lorentz force 
and the magnetic energy density as anti-correlated [8J, J9| , 
however, the Lorentz force and the magnetic energy den- 
sity correlation can be calculated, as for the their auto- 
correlations, as pointed out in Ref. [lOl- The cross- 
correlator between pB and Lb is given by: 

{pBik)LBik'))^^-^^JdpPB{p)PBi\k-p\) 

X (1- 1(72 +/32)+ 27/3^ -^2)(7) 

The above formula can be computed analitycally and in 
the appendix we show the analytical results both for 
Ub — —2.5 and ns — 2. The general behavior in the 
large scale limit of the spectrum (for k « ko), the 
range relevant for CMB anisotropics, is: 

{pB{k)LB{k)) = -i|ps(fc)|' forns > -3/2 

{pB{k)LB{k)) = -C|pB(fc)|2 forns < -3/2 

where C ^ 0{1)- In the appendix we show that the exact 
evaluation of the cross correlation has a small effect - in 



particular for strongly negative ub - with respect to our 
choice of considering the Lorentz force and the energy 
density as fully anti-correlated. 



B. Magnetized Scalar Initial Conditions with 
matter Corrections 



The second issue is related with the initial conditions 
for magnetized scalar perturbations. In our previous 
works [Sl, |9[ we computed the intial conditions in the ap- 
proximation of a universe dominated by radiation, follow- 
ing the results presented in Ref. 130| we investigated the 
importance of matter contributions to initial conditions. 
We computed the initial conditions with the inclusion of 
matter corrections and we give the details in appendix 
C. The results show that the introduction of matter cor- 
rections in the initial conditions have a negligible impact 
on the numerical computation of the magnetized scalar 
CMB anisotropics. 



III. CONSTRAINTS FROM CURRENT CMB 
DATA 

We derived the constraints on PMF with current CMB 
anisotropy data performing an analysis of the WMAP 
7 year [2l, [2^^.CBAR [22] data in temperature and 
of the BICEP [H and QUaD ^ data in polarization. 
In order to decrease the correlations between different 
data sets which cover the same region of the sky, we 
proceed as in Ref. [3l|. We remove in the analysis the 
following CMB band powers: a) all the QUaD TT band 
powers since they overlap with data from the 'CMB8' 
region of ACBAR, b) the ACBAR band powers with £ < 
910 and i > 1950 to avoid overlap with WMAP (which 
is cosmic variance limited up to £ — 919 [25|, l26[) and 
contamination from foreground residuals, respectively, c) 
the QUAD TE band powers which overlap with WMAP 
ones, the QUAD EE band powers which overlap with 
BICEP, d) the BICEP TT, TE band powers (i.e., we use 
just EE and BB information from BICEP). 

We use a modified version of CosmoMC 20] in order to 
compute the Bayesian probability distribution of cosmo- 
logical parameters, including the magnetic ones. We vary 
the baryon density Wf, = flbh'^, the cold dark matter den- 
sity ojc = ric/i^ (with h being ifo/100kms~^Mpc~"^), the 
reionisation optical depth r (not to be confused with the 
conformal time r), the ratio of the sound horizon to the 
angular diameter distance at decoupling 9, \n{10^'^ As), 
ns and the magnetic parameters -B^Mpc (in units of 
lOnG) and ub- As priors we use [0 , 10] for _BiMpc/(10nG) 
and [—2.9,3] for ub- The damping scale kjj defined in 
Eq. ([3]) is obtained as a derived parameter (in units of 
Mpc^^), as well as Hq. 

We assume a flat universe, a CMB temperature 
Tqmb = 2.725 K and we set the primordial Helium frac- 
tion to 2/Ho = 0.24. We use the recombination routine 
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TABLE I. Mean parameter values and bounds of the central 
95%-credible intervals without (left column) and with (right 
column) PMF. Note that Hq is a derived parameter. 



Instrument 


LFI 


HFI 1 


Center frequency GHz 


70 


100 


143 


Mean FWHM (arcmin) 


13.0 


9.6 


7.0 


AT/r per FWHM" (Stokes /) 


4.45 


2.12 


1.56 


AT/T per FWHM" (Stokes QSzU) 


6.29 


3.39 


2.90 



TABLE II. Planck channels considered in this analysis and 
relative performances in the extended mission (four surveys). 



RECFAST vl.5 [SJ]- We restrict our analysis to three 
massless neutrinos (a non- vanishing neutrino mass leads 
to a large scale enhancement in the power spectrum of 
CMB anisotropics in the presence of PMF |lfl] and would 
not change our results). The pivot scale of the primor- 
dial scalar was set to fc, = 0.05 Mpc^^. In order to 
fit WMAP 7 years, ACBAR and QUaD data, we use 
the CMB power spectra and we follow the method im- 
plemented in CosmoMC consisting in varying a nuisance 
parameter Agz which accounts for the unknown ampli- 
tude of the thermal Sunyaev-Zeldovich (SZ) contribution 
to the small-scale CMB data points assuming the model 
of [37|]. We use a CAMB accuracy setting of 1. We sam- 
ple the posterior using the Metropolis-Hastings algorithm 
[38| at a temperature T = 2 (the temperature parameter 
in the CosmoMC code is used to sample the probability 
P as P^/^, allowing a better exploration of the distribu- 
tion tails), generating four parallel chains and imposirig 
a conservative Gelman-Rubin convergence criterion ^0^ 
ofR-1 <0.01. 
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UJb 


0.0227 
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UJc 
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0.089 
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log [lO^^^s] 


3.1 


3.08«:g? 


ns 


0.960 


0.961 ± 0.008 


-BiMpc/nG 




< 2.7 


TIB 




< -0.05 


i?o/kms-^Mpc-^ 


72.4 


72.3 ±1.3 



TABLE HI. Input and mean parameter values with bounds 
of the central 95%-credible intervals for the Planck simulated 
data. Note that Ho is a derived parameter. 



In Table U are reported the results of our analysis on 
current CMB data. We compare the results for the six 
standard cosmological parameters both considering or 
not the PMF contribution. The comparison shows that 
neither the means nor the bounds of the cosmological 
parameters of the ACDM model with reionization are 
basically affected by the presence of PMF: this means 
that Bx and ns are not degenerate with the other six 
parameter of the concordance cosmological model. 

In Figs. S and [5] we show the plots of the MCMC 
results. We derive the following constraints on the am- 
plitude and spectral index of PMF: i?iMpc < 5.0 uG and 
n < —0.12 at 95% confidence level. 

In Fig. \S\ we compare the results obtained includ- 
ing or excluding the scalar contribution. We note how, 
as expected from previous results on CMB anisotropics 
@> fi ) the scalar contribution has a very little impact 
on the constraints on the amplitude Bx and spectral in- 
dex ub ■ This is due to the different shape of the CMB 
temperature spectra induced by the scalar and vector 
contribution of a stochastic background of PMF: com- 
pared to the CMB anisotropics sourced by the adiabatic 
mode, the PMF scalar mode is the dominant contribu- 
tion on large scales, whereas the vector is the important 
one at high multipoles. Since the PMF contribution does 
not suffer of the Silk damping, only the vector contribu- 
tion has the main chance to affect the CMB temperature 
sourced by the standard adiabatic mode. 

A stochastic background of PMF has also an impact 
on the matter power spectrum as shown in [8|, |lO|. In 
order to investigate if the addition of matter power 
power spectrum data could improve the constraints on 
the PMF, we perfomed a MCMC analysis adding the 
large scale structure (LSS) data of SDSS LRG DR4 
|34| to the CMB anisotropy data used before. We 
used the DR4 release because the DR7 requires the 
use of non-linear tools to compute the matter power 
spectrum which are unavailable for ACDM plus PMF. 
We have verified that the addition of the SDSS data 
does not introduce any significant improvement on the 
constraints on PMF amplitude. The reason for the 
lack of improvement is that the matter power spectrum 
data reach scales of the order of 0.2 hMpc~^ where only 
PMF with positive spectral index have an impact, but 
PMF with blue spectral indices are already strongly 
constrained with CMB anisotropy data, therefore the 
results do not show any significant improvement. This 
result of a negligible impact of SDSS DR4 data with 
respect to current CMB data is in agreement with |23| . 

Let us compare our constraints with other investiga- 
tions. After our manuscript appeared on the archive, 
Shaw and Lewis obtained similar CMB constraints - of 
the order of 6 nG at 95% confidence level - with simi- 
lar data sets in Ref. [231. ^^ UM stronger constraints of 
the order of 3 nG at 95% confidence level were reported. 
The reason for slightly different result might be either 
in the different datasets or in the different methodology 



used. Concerning the use of matter power spectrum, the 
present manuscript and Ref. [23l | agree that LSS data 
add very little to CMB current constraints, when only 
linear scales are correctly included: we obtain analogous 
resuhs either with SDSS DR4 i3j] or 2dF [35.], without 
changing the CMB constraints in Table I. Another differ- 
ence might be in the use of CBI data [33| in [24I , whose 
foreground contamination at high £ has been stressed in 
Ref. [231 : we believe that the difference might be due to 
either methodology and/or to the different use of CMB 
data. 

Let us end this section by commenting on the con- 
straints on the magnetic spectral index we derive, i.e. 
riB < —0.12 at 95% confidence level. This would im- 
ply that positive values for ub are disfavoured and so 
would be causal mechanism which would have produced 
such PMF [42] . Certainly there is a trend in Bx , ns 
which allows for larger B\ for negative values of ub'- if 
we fix riB to the causal value of 2 we obtain a much 
stronger constraint on the amplitude, i.e. SiMpc < 0.023 
nG 95% confidence level. At the same time, we have 
verified that for undectatable values by CMB of Biupc - 
fixed to 10~^ nG - no constraints are derived for ns, as 
expected. The above results are obtained by the assump- 
tion of a linear prior on -BiMpc- We have also performed 
just for comparison cosmological parameter extraction by 
sampling Logiol-BiMpc/(10nG)] in the range [—7,1] [4j] 
and we show the posteriors in Fig. [6l Since the log 
prior overweights the importance of small amplitudes for 
^iMpc compared to the linear prior, the constraints ob- 
tained on SiMpc by sampling logarithmically are tigher 
- 1 nG at 95% confidence level - , but depend strongly 
on the interval. As byproduct of the importance given to 
small amplitude for J3iMpc in sampling logarithmically, 
Ub is unconstrained. Similar prior issues occur for other 
important cosmological parameters which are under the 
threshold of detection, as the amplitude of g ravitational 
waves in standard inflationary scenarios 44 [ . 



IV. FORECASTS FOR PLANCK 

In the perspective of the forthcoming Planck [2J] 
data we performed a MCMC analysis using simulated 
Planck data. We created the mock data considering 
the updated angular resolution and sensitivities [40, |4l| 
for the extended and approved mission duration of four 
sky surveys, which are summarized in Table HIl We used 
the combination of the three central frequency channels 
where the CMB is dominant with respect to foreground 
contributions: 70 GHz, 100 GHz, 143 GHz, which have 
been combined with the inverse noise variance weighting 
technique. We used the 92 % of the sky and assumed 
only CMB neglecting all possible foreground contamina- 
tion. In Table IIIII we report the results for cosmologi- 
cal parameters with the input parameters of the fiducial 
model which we choose very close to the WMAP 7 yrs 
best fit model without PMF. 



In Fig. [7] we show the comparison between the pa- 
rameters constrained by current CMB data and the ones 
which win be constrained by Planck: note the great im- 
provement given by Planck on the constraints on cos- 
mological parameters in presence of PMF. We forecast 
the fohowing constraints by Planck: SiMpc < 2.7 nG 
and ub < —0.054 at 95% confidence leveL The improve- 
ment obtained by Planck alone on the constraint on 
Bx is mainly due to the better measurement of the CMB 
temperature anisotropies at high multipoles. The factor 
of 2 improvement in the constraints on B\ reached by 
Planck corresponds to a factor 16 improvement in the 
characterization of Cf. 



the initial conditions for magnetic scalar mode in pres- 
ence of matter corrections showing how these corrections 
do not affect the results contrary to what claimed in [30| ■ 
On the basis of previous works, we have considered 
only the scalar and vector contribution, and by using 
their regular initial condition, we constrain i?iMpc < 5.0 
nG and ub < —0.12 at 95% confidence level, with 
the most updated combination of CMB anisotropies. 
Planck will be able to constrain the spectrum of a 
stochastic background of PMF even further at the level 
of 2.7nG. 



V. CONCLUSIONS 

We studied the constraints on a stochastic background 
of PMF by current and forthcoming CMB data. In doing 
this, we have improved the theoretical CMB predictions 
in presence of a stochastic background of PMF. We gave 
approximations for the relevant components of the EMT 
for any spectral index ub ■ We considered the correlation 
between ps and Lb and showed how the previous choice 
of total anti-correlation [^, Q was a rather good approxi- 
mation, in particular for red values of nB ■ We computed 

I 
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Appendix A: PMF EMT spectral fits 

Since the spectral behaviour is rather complicated, in order to have the best fit possible we decided to divide the 
fits for spectral index ranges. Together with the natural separation between indices smaller and greater than —3/2, 
which is necessary due to the completely different behaviour, we added a further separation between negative and 
positive spectral indices. The leading terms remain the same as in the infrared limit of the exact spectra: white noise 
for blue indices and infrared dominated for the red ones. 
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1. Scalar Spectra 
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FIG. 1. Comparison of spectral fit (dashed line) and exact spectrum (solid line) for the rescaled magnetic energy den- 
sity |ps(fc)p/(7 (top panel), scalar Lorentz force |Ls(A:)p/t/ (middle panel) and the vector part of the anisotropic stress 
|n<^'(A:)|V(2(7) (bottom panel) as a function of k/ko with U = (A^fc^"^+*)/(5127r'*). In the left (right) column hb = 2.3 
(hb ~ —2.7) is displayed. 



Coefficients 

For positive spectral indexes: 
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for the negative spectral indexes we have: 
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2. Lorentz Force Spectra 
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Coefficients 

For positive spectral indexes: 



Af ==0.933635 



0.00460612 

riB 



0.505278nB + 0.183487n?, - 0.0238037n?, - 0.00985191ni + 0.00437658ni - 0.000504247n 



A^ = 0.22309 - ^'°^^^^^ - 0.152155nB + 0.427087n| - 0.184484n| - 0.0111374n| + 0.0292611n| - 0.00571069n' 



A^ = 1.84015 



riB 

0.319013 9 c, 4 . fi 

3.60452nB + 2.88574no - 0.797507n?j - 0.145007ni + 0.116527n| - 0.0163659n| 

riB 



for the negative spectral indexes we have: 
1 



Bt 



B^ 



41 



(1630 + 4240nB + 3360n| - 2080n|j - 1960ni + 1970n| + 1559n|) 
100 25n_B 

-—(-854 - 2838nB - 2710n| + 1390n?j + 1705nl - 1530n| - 1340n?j) 

100 oub 



for strongly negative: 
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3. Vector Spectra 
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Appendix B: Cross Correlators Exact Solutions 



To give the reader an idea of the complete analytical form of the Lorentz force and energy density cross correlator 
we show the exact results for two representative spectral indices ub = 2 and ub — —2.5: 



{pB{k)LB{k))\nB=2 
(pB(fc)ii3(fc))|„B = -5/2 



10247r5a8 
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105fc3 
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where k — k/kjj. In Fig. [2] we show the comparison between the energy density the Lorentz force and their cross 
correlation for ub — —2.5 and ub = 2: the cross-correlation is negative in the whole range of scales. In Fig. [3] we 
show the impact of various assumptions for the cross-correlation between the magnetic energy density and the Lorentz 
force on the power spectrum of CMB temperature anisotropics. 




2.0 0.0 



FIG. 2. The cross-correlator k^ {pB{k)LB{k)) /U (dotted line) for ub ~ 2 and ub = —5/2 is plotted versus k/ko in comparison 
with fc^ipfl(fc)|V[/ (solid line) and k^\LB{k)\^ /U (dashed line) with U = (A2fc^"^+^)/(5127r'*). 
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FIG. 3. We computed the scalar power spectrum with the contribution of the Lorentz force-energy density cross correlation for 
Ub ~ 2 (top panel) and ns = —2.5 (bottom panel). Solid line represents the correct assumption on the cross-correlation, the 
dotted line represents the result considering the Lorentz force and the energy density fully anti-correlated, triple dotted-dashed 
line represents the uncorrelated sum and the dashed line represents the result assuming full correlation. 



Appendix C: Scalar Initial Conditions with Matter Corrections 



We computed the initial conditions for magnetized scalar perturbations for the leading regular growing mode which 
is the one which contributes to CMB anisotropics, we do not consider any isocurvature or decaying modes. Initial 
conditions are the solution of the coupled Einstein Boltzmann equation system with PMF contributions on long 
wavelenghts at early epochs. The Einstein equations with the contribution of PMF in the synchronous gauge (with 
the notation of [3a]) are: 

1 



f^^V - ^'Hh ^ -47rG'a^(I]„ p„(5„ -I- ps) , 



h + 2nh-2k^r] = -SttGc^ I]„ cj „p„(5„ 



3 / 



h + 6ij + 2n{h + 6?)) - 2k'^ri = -247rGa^[S„(p„ + P„)^n + ctb], 



(CI) 



where n represents the various species of the plasma: b for baryons, c for cold dark matter (CDM), 7 for photons 
and I' for massless neutrinos. The scalar metric perturbations are represented by the two scalar potential h, rj while 
fluid perturbations are represented by Si, 9i, Gi which are respectively the density contrast the fluid velocity and the 



10 



anisotropic stress in the notation of (36|. The overdots denote the derivative with respect to conformal time. On 
long wavelenght we can expand the metric and fluid perturbations in series oi kr << 1. In our previous works [8|, |9[ 
we computed the intial conditions deep in the radiation era with the approximation of a universe dominated only 
by relativistic degrees of freedom (radiation and massless neutrinos). Usually considering a{T) ex r is a rather good 
approximation for setting the initial conditions (see however Ref. ^30] for a different claim stressing the importance of 
matter corrections) . We have therefore extended our previous results [9| to the case in which the matter contribution 
is taken into account. Considering the matter contribution in the initial condition needs the introduction of matter 
correction also in the evolution of the scale factor 451: 
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flft and the Hubble parameter can be expanded as H oc ^ 
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J. The Einstein equations become: 

(C2) 
(C3) 
(C4) 



■^Pu^v + {pv + p^)y— + LB 



(C5) 



where PMF contributions are given in terms of their ratio with the fluid radiation density: VLb — ps/ Prad^ Lb = 
Lb / Prad and gb = obI Prad (with prad = Pu + Py) ■ Thc initial conditions at leading order with PMFs which include 
the matter corrections are: 
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where R^, = Pvl{pi, + p-y). We note that our results are in agreement, within our notation, with the one presented in 
[lO| . We verified that the inclusion of the matter correction in the initial conditions in our modified CAMB code does 



not produce any appreaciable change in the results: /S.C1/C1 
with Ref. ;30]. 
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10 ^ which is the level of numerical noise, in contrast 
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FIG. 4. Results of the MCMC constrained with WMAP 7 years, BICEP, QUaD and ACBAR data. Curves are the 68% and 
95% confidence leveL Note that B\ (with A = 1 Mpc) is in 10 nG units. 
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FIG. 5. Results of the MCMC constrained with WMAP 7 years, BICEP, QUaD and ACBAR data. Curves are the 68% and 
95% confidence leveL Solid lines are the results considering both scalar and vector magnetic modes while dashed lines are the 
results taking into account only the vector contribution. Note that Bx (with A = 1 Mpc) is in 10 nC units. 
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FIG. 6. Results of the MCMC for the PMF spectral index constrained with WMAP 7 years, BICEP, QUaD and ACBAR data 
and logarithmic prior for the PMF amplitude. Curves are the 95% and 68% confidence level. 
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FIG. 7. Comparison of the results of the MCMC with real data used in Figs. 5, 6 (solid line) and simulated Planck data 
(dashed line). Curves are the 95% and 68% confidence level. Note that Bx (with A = 1 Mpc) is in 10 nG units. Note that ko 
is a derived parameter. 
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